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In this manuscript the B —>■ 'ktVt decay is investigated in the context of the 2HDM type II 
extension of the SM. In particular, a prediction for the ratio of semileptonic branching fractions 
from f = r and i = e, fi is produced and an exclusion based on the MSSM parameters tan /3 and 
from a recent Belle measurement is reported. An alternative variable to the ratio of branching 
fractions that could offer increased sensitivity is discussed that could result in increased sensitivity 
to rule out contributions from additional scalar mediators. 


I. INTRODUCTION 

In recent measurements the i3-Factory experiments 
BaBar, Belle and LHCb have reported large disagree¬ 
ments in semi-tauonic decays involving ratios of charmed 
final states m ■ A recent average from Ref. [5] quan¬ 
tifies the overall disagreement with the Standard Model 
(SM) expectation at the 3.9 a level. This disagreement is 
an interesting anomaly that deserves more experimental 
and theoretical study, in particular due to the absence 
of any clear new physics signal from the experiments 
at the energy frontier of particle physics. An interest¬ 
ing alternative probe to see similar effects is given with 
B -A TTTVr'- although CKM suppressed, the B ^ ttt Dr 
decay is less phase space suppressed with respect to the 
transition involving light leptons £ = e, p, as its charmed 
partner decays, resulting in a high fraction oi B -A ttt Dr 
over B > 7 r£i/f in the SM. In this manuscript a predic¬ 
tion for the ratio 


struction of the four-momentum transfer squared and the 
manuscript concludes in Section |VI| with a summary of 
the key results. 

II. THE B ^ -IT r Dr DECAY IN THE SM 

The effective SM Lagrangian describing the b ^ u£ D^ 
transition is well known in the literature and given by 

Aff = Kf, {uinPLb) (D^PLi) + h.c. (2) 

with the projection operator Pj, = (I — 75 ) /2 and Gp 
Fermi’s constant. The B -A iriDi decay amplitude de¬ 
pends on one non-perturbative hadronic matrix element 
that can be expressed using Lorentz invariance and the 
equation of motion in terms of P —>■ tt form factors. The 
P —> TT transition form factors /+/o are defined by 


Rr = 


B{B^ 


B{B 


( 1 ) 


for the SM is given and its modifications in the con¬ 
text of the two-Higgs doublet model (2HDM) of type 
II are discussed. A first exclusion of the MSSM param¬ 
eters tan/3 and is carried out using recent limits 

on the P -A TTT Dr branching fraction reported by the 
Belle experiment [S]. As will be shown. Equation can 
be predicted with a precision of a few percent as many of 
the theoretical uncertainties associated with the hadronic 
form factors cancel and no dependence on the CKM ma¬ 
trix element \Vub\ remains. Similar predictions have been 
discussed in Refs. [iHn]. Recent progress in the under¬ 
standing of the P —>■ TT form factor as well as the availabil¬ 
ity of first experimental results give a strong motivation 
to revisit this topic. 

This manuscript is organized as follows: Section |ll] 
briefly reviews the P -A ttt Dr decay in the SM and sum¬ 
marizes the state-of-the-art knowledge of the P —)■ tt form 
factor. Section im discusses the modifications of the de¬ 
cay rate in the context of the two Higgs dou blet model 
(2HDM) type II interactions. In Section IV the predic¬ 
tion of the ratio of branching fractions from B ^ tt t Dr 
and P —iriDg is discussed and the exclusion of the 
2HDM type II parameter space is presented. Section [V| 
discusses a more sensitive variable relying on the recon- 
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with q = p — pr denoting the four-momentum trans¬ 
fer between the P-meson and the final-state pion of the 
semileptonic decay. Further, rriB denotes the P-meson 
mass and rrir the pion mass. The differential decay rate 
as a function of with its full lepton mass dependence 
is given by 


dr(P Tr£Ui) _ 8 \pr\ G]r \Vub\ 
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with the Helicity amplitudes Pq and Ht and \pr \ the ab¬ 
solute three-momentum of the final state pion. The ab¬ 
solute three-momentum is related to the four-momentum 


^ In this section i denotes any massive lepton unless stated other¬ 
wise. 













2 


transfer squared as 
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Setting the lepton mass mi to zero one recovers the 
expression 


AT{B ^ t: e.vi) _S\p^\G], \Vub\ q 
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that holds to be an excellent approximation for mi = m^ 
or mg = TO^. The helicity amplitudes ifo/t are related to 
the form factors defined in Equation as 



(7) 


(8) 


The form factors /+/o are functions of the four- 
momentum transfer squared and need to be calculated 
with non-perturbative methods. 

Analyticity and unitarity impose strong constraints on 
heavy meson decay form factors [T2HI6] . In the following 
the form factors are parametrized using a series expan¬ 
sion what involves a mapping of the variable to the 
variable 


z{t,to) 


y/t+-q^ - y/t+ - to 
\/t+ - g 2 + ^t+ - to ’ 


(9) 


with 


t± = {ms ± , to = (ms -I- m^) (v^ms - y/m^f ■ 


In this expansion the form factors /+ and /o are assumed 
to be analytic in z expect for a branch cut in [t+, oo) and 
poles in [t _, t_|_]. Using Blaschke factors these poles can 
be removed, c.f. Ref. [HIIE], and the form factors can 
be expanded at[^ 
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with TOg = 5.325 GeV the i3*-meson mass, and the 
number of expansion parameters Equation 
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IS 


2 This expression is recovered by starting with an expansion 
P{z)4>{z) f = , with Pi the Blaschke factor and by 

setting the outer functions (p to unity. 


Lattice bi bl b\ bi 


0.407(15) 

-0.65(16) 

-0.5(9) 

0.4(1.3) 

/i=0 

0.507(22) 

-1.77(18) 

1.3(8) 

4(1) 

Lattice-!-Exp 

bh 

b{ 

h\ 

bl 


0.419(13) 

-0.495(54) 

-0.43(13) 

0.22(31) 

/i=0 

0.510(19) 

-1.700(82) 

1.53(19) 

4.52(83) 


TABLE 1. The BCL expansion parameters of m are sum¬ 
marized: ’Lattice’ shows the lattice only results and ’Lat- 
tice-|-Exp’ the result of simultaneous fit the lattice informa¬ 
tion and all suitable measured partial branching fractions of 
B —>■ nii/i with i = e,fi. The full correlation information 
between the parameters can be found in Tables XIV and XIX 
of Ref. [IT] . 


known as the Bourrely-Caprini-Lellouch (BCL) expan¬ 
sion |16] in the literature. 

Ref. [ITj recently reported (2-|-l)-flavor lattice QCD 
calculations of the two form factors /+ and /o in the 
BCL expansion. In addition. Ref. HU also carried out a 
global analysis of the lattice QCD expansion parameters 
and all suitable B —> tt iDg partial branching fractions 
with the goal to extract \Vub\- The post-fit expansion pa¬ 
rameters from the combined analysis represent the cur¬ 
rent best knowledge of the B ^ tv form factor, under 
the assumption that no new physics is affecting the light 
lepton decays. Both sets of expansion parameters are 
summarized in Table |T] and Eigure[^ (left) shows the pre¬ 
dicted differential rates for R —>■ tt r and B —b iriDi 
using the expansion parameters from the global analysis. 


III. MODIFICATIONS IN THE 2HDM TYPE II 

In the two-Higgs doublet model extension of the SM, 
which describes the Higgs sector of the Minimal Super- 
symmetric model at tree level, the B —>■ tt t Dt decay 
also can be mediated by a Higgs boson. Such an 
additional scalar mediator has drastic consequences as 
due to the Spin I nature of the boson the i? —)• tt 
transition occurs predominantly as a P-wave and is ac¬ 
cordingly suppressed by a factor near the maximal 
four-momentum transfer squared (where |p(r| ~ 0). The 
B ^ TT form factors reach their maximal value at the 
maximal four-momentum transfer and non-P-wave con¬ 
tributions from e.g. a scalar Higgs boson would greatly 
affect the rate in this part of phase space. 

The contributions of the charged Higgs Boson to P —>■ 
TTTVr decays can be incorporated into Equation]^ by the 
replacement Bin] 

uSM , u2HDM uSM f ^ P 

‘ ‘ ‘ 1 ml^ l-mjmb 

( 12 ) 

where tan /3 is the ratio of the vacuum expectation values 
of the two Higgs doublets and is the mass of the 
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FIG. 1. (left) The differential rate for B ^ n ivt with £ = e,fi and B —>■ tt r i/r in the SM with associated errors are shown. 
As input for the form factors the ’Lattice+Exp’ values of Table [l| are used (right) The differential rate of B —>■ tt r i/r is shown 
for the SM and two MSSM working points: tan P/mff+ = 0.3 (B'(^°^ < Hf^) and tan/3/m^+ = 0.4 > Hf^). 


charged Higgs boson and Hf^ = Ht , cf. Eq. Further 
mulrrif, is the ratio of the u- and &-quark masses at an 
arbitrary (but common) mass scale. This factor leads 
to a negligible modification and can in principle safely 
be ignored. Equation proofs sufficiently accurate to 
probe charged Higgs bosons with a mass of more than 
about 15 GeV, what is the experimental interesting re¬ 
gion as mass scales below 15 GeV are excluded by 5 —>■ s 7 
measurements already [ 20 ] or other constraints 

Figure]^ (right) compares the SM differential rate with 
two working points for tan /3 and mH+. The rate is 
only sensitive to their respective ratio and for a value 
of tan P/mH+ = 0.3 the negative interference causes 
B( 2 HDM smaller than Hf^, resulting in a lower 

B — )■ tttDt rate than for the SM. For tan l3/mH+ = 0.4 
the additional scalar transition starts to dominate the 
matrix element, modifying the differential decay rate 
strongly near the kinematic endpoint. For larger values 
of tan f]/mH+ the 2HDM contribution completely domi¬ 
nates the total and differential rate. 


IV. PREDICTION FOR AND EXCLUSION 
FOR tan /3 AND 

By integrating the differential rates Equations and 
over the allowed kinematic range in q^, a prediction for 
Rtt can be obtained: 

^ _ r(B^7rTF^) _ dg2 dr(B ^ ttt t>^)/dg2 

r(i? —>■ nivi) dq'^ dr(i? —)■ n£Pi)ldq‘^ 

(13) 

with ~ and the light lepton 

masses can be set to zero. The CKM matrix element 
\Vub\: Fermi’s constant, as well as other constant terms 
cancel in the ratio. Using the B —>■ tt BCL parameters 


from the global analysis of Table |T] the following predic¬ 
tion for the SM value is obtained: 

= 0.641 ± 0.016 . (14) 

The rates in the numerator and denominator are strongly 
correlated (correlation p = 0.94) if the integration is car¬ 
ried out over the full kinematic range, resulting in the 
large desired cancelation of non-perturbative uncertain¬ 
ties from the form factors. This number can be com¬ 
pared to the recently reported first measurement from 
Belle [5]0 

^meas ^ 105^0.51. (15) 

Although statistically limited, this ratio already has some 
power to rule out several working points in tan /3 and 
m//+ as many predict a very strong enhancement of 
This can be seen in Figure (left) from the predicted 
(dark grey) and observed (light grey) B,r values. 

To perform a proper exclusion in terms of tan/3 and 
mH+ though the impact on the acceptance and selec¬ 
tion efficiencies should be studied. This can be done for 
instance by reweighing the used SM Monte Garlo (MC) 
samples to the probed MSSM working point with weights 
such as 

_ Fmc dF(tan/3,77i.g+)/dg2 

^ r(tan/3,mjj+) dFMc/dg^ 

where the first factor ensures that the branching frac¬ 
tion in the simulation is unaltered and the quantities are 
Fmc the total rate in MG, r(tan/3,m//+) the total rate 


® The reported value for the branching fraction of Ref. [6] was di¬ 
vided by the current world average of the B® ui branch¬ 

ing fraction of (1.45 ± 0.05) X 10“"^ from Ref. 1231 . 
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FIG. 2. (left) The evolution of R-^ as a function of tan/3/m£f+ is shown (dark grey). The experimental 68% Cl of Ref. [6] is 
shown (in light grey) and MSSM parameter points that predict ratios larger than ~ 0.48 can be excluded at 95 CL under the 
assumption that efficiencies and acceptances are modified negligibly with respect to the current statistical precision, (right) 
The 68% (light grey) and 95% (grey) excluded MSSM parameter points in the tan j3 — mu+ plane. 


at a given working point in tan /3 —171 ^+, and dF/ dq^ the 
corresponding differential rates. In addition the fraction 
of negative and positive helicity r-leptons has to be cor¬ 
rected for as e.g. subsequent leptons are either emitted 
preferentially in or opposite to the t direction. The evo¬ 
lution of the helicity fraction as a function of tan/3/m//+ 
is briefly discussed in appendix 

Assuming that such effects are small with respect to 
the current statistical precision, a first exclusion can be 
carried out: Figure (right) shows the excluded region 
of tan/3 — mH+ at 95% CL (dark grey) and 68 % CL 
(light grey). A comparison with the current bounds from 
B ^ T Dr can be found in appendix 


with a lower cut on q^. This is motivated that in 
the 2HDM type II for large values of ta.n P/mH+ the 
constructive interference is dominating the rate Equa¬ 
tion]^ and due to the absence of the P-Wave suppression 
strongly modifies the large q^ region. A fully differential 
analysis of course contains the most power to probe and 
distinguish such a scenario with a scalar mediator from 
the SM, but with the current small experimental sensi¬ 
tivity a measurement that probes the high q^ range in 
one bin already could lead to an improved rejection of 
large tan P/mH+ ratios. For instance: the SM prediction 
for a measurement in [(?max/ 2 , 9 max] 


V. OTHER OBSERVABLES WITH INCREASED 
SENSITIVITY 

An interesting alternative observable besides Rr could 
be the ratio of partially integrated and measured rates 


dSM / 2 _ 2 

■^TT VTmin ^max 


/ 2 ) 


clF(i3 ^TTT Dr)/dq^ 

^min ^max/ ^ _ 

/2 d?" dTiB^ire D,)/dq^ 

'^min 'imax/ ^ 


1.012 ±0.008. 


(17) 


and the ratio of branching fractions at for example 
tan l3/mff+ = 0.4 results in a prediction about a factor 
of two larger for Rr{q^^^/2) = 2.09 ± 0.03. In compari¬ 
son, the fully inclusive rate ratio for the same tan/3/m//+ 
working point results in R^ = 1.01 ± 0.04 and shows 
a smaller enhancement with respect to the SM value 
Eq. |14[ The reduction in experimental sensitivity on R^ 


due to only analyzing half of the allowed phase space of 
course also has to be taken into account: a rough estimate 
for the increase in statistical error, assuming that back¬ 
ground is distributed uniformly in q^ and that the signal 
is SM like, is a factor of about « y/2. These numbers 
imply an overall improved sensitivity for a measurement 
in [(Imax/^i 9 max] ^^e fully inclusive R^ value to ex- 
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elude tan/3/m//+ working points. In addition, as pointed 
out already by Refs. [3 El], the cancellation of the non- 
perturbative error further improves allowing for a more 
precise predictions of the ratio of partial rates. 


VI. SUMMARY AND CONCLUSION 

The B ^ TT T Vr decay offers an interesting alternative 
decay to probe the deviations observed in semi-tauonic 
decays with charmed final states. In this manuscript a 
prediction for the ratio of semi-tauonic and light lepton 
total rates, is presented using the state-of-the-art 
knowledge of the B —>■ tt form factor. In the SM a value 
of i? 7 r = 0.641 ± 0.016 is found. The impact of the pres¬ 
ence of charged Higgs boson contributions in the context 


of the 2HDM type II is discussed and a first preliminary 
exclusion in tan /3 — m}j+ using a recent Belle measure¬ 
ment is carried out. In addition the idea is presented to 
measure in future measurements the high g^-region, as it 
has a higher sensitivity to probe the MSSM parameter 
space in the 2HDM type II model. 
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Appendix A: Helicity fractions in the 2HDM type II 

In the SM the t~ polarization is about 35% positive 
and 65% negative. This relative fraction changes dramat¬ 
ically in the 2HDM type II: in the range of tan j3/mH+ € 
[0,0.25] the negative helicity fraction increases to almost 
85%. This has drastic consequences if the r“-leptons are 
reconstructed using e.g leptonic channels, t~ —>■ t~ vf. 
the secondary lepton spectrum in the B-meson rest frame 
will become harder with respect to the SM spectrum as 
the lepton is preferably emitted in the direction of the 
T“-lepton altering the acceptance and selection efficien¬ 
cies. After this turning point the positive helicity starts 
to dominate and at tan/3/m//+ = 0.7 about 95% of the 
r“-leptons have positive helicity. If reconstructed again 
via T~ —>■ £~ Vt V£, the secondary lepton spectrum now 
is softer than the SM spectrum as the leptons preferably 
are emitted opposite to the r”-lepton direction. Fig¬ 
ure (left) shows the helicity fractions as a function of 


tan /3/ m}j+ for B —>■ tt r Ur- 


Appendix B: Comparison with B ^ tVt 

Figure (right) compares the 95% exclusion limits of 
the 2HDM type II parameter space from B —ttt D^. 
with the exclusion from B — >■ t Dr'- Both decays involve 
the same quark lines and the existing measurements for 
B ^ T Dr result in tighter constraints as the individual 
B —> TT T Dr branching fraction. The two red regions is 
from B ^ T Dr and the dark red region (dashed black 
boundary) is from B ^ ttt Dr- As input a world aver¬ 
age of B{B -D- rDr) = (1.03 ±0.20) x 10“^ was used, 
obtained from averaging the latest Belle results [l5l 
with the values reported from BaBar m- Further, for 
\Vub\ and the B-meson decay constant the reported val¬ 
ues of Ref. m and [28] were used, respectively. 


r Helicity Fraction 
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FIG. 3. (left) The T“-lepton positive (+) and negative (-) helicity fractions as a function of ta.n jS/m^+ for B — >■ 'kt~ Vt 
is shown, (right) The 95% exclusion limit from Ref. [6] is compared with the exclusion limits from B ^ t Vr- As input an 
averaged branching fraction of B{B — > r hr) = (1.03 ± 0.20) x 10“"* and \Vub\ = (3.72 ± 0.16) x 10“^ was used. 






